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SUMMARY 


The  principle  uncertainty  in  the  visco-plastic  model  of  hypervelocity  impact 
is  the  lack  of  data  on  strain-rate  effects  under  the  severe  conditions  of  interest. 
Estimates  of  /x  for  iron  and  aluminum  are  first  determined  that  are  compatible 
both  with  observed  strain-rate  effects  in  plastic  waves  and  the  very  abrupt  shock 
profiles  known  to  exist  at  pressures  greater  than  about  .2  megabar. 

In  the  second  part  of  the  report  detailed  calculations  are  presented  describing 
the  cratering  process  for  a  cylindrical  projectile  impacting  a  thick  target  of  like 
metal.  It  is  shown  that  momentum  scaling  cannot  be  justified;  energy  scaling  for 
geometrically  similar  impact  situations  is  accurately  predicted  upon  taking  into 
account  the  different  rates  of  shock  propagation  during  the  early  stages  of  the 
cratering  process.  The  depth  of  penetration  is  therefore  related  to  impact  velo¬ 
city  Vq  according  to  Pg  =  KgV^  2/3  for  Vq  =  Vq*.  Strength  effects  of  the  target  do 
not  vanish  with  increasing  impact  velocity. 

In  the  third  part  of  the  report  the  cratering  and  the  penetration  process  for 
thin  plates  is  studied.  PICWICK  calculations  are  presented  and  the  results  dis¬ 
cussed  in  terms  of  the  interactions  of  shock  waves,  rarefaction  waves,  and  con¬ 
figuration  geometry. 
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INTRODUCTION 


Previous  reports  have  described  the  visco-plastic  model  that  has  evolved  to 
describe  the  hypervelocity  impact  cratering  and  penetration  process  [1,  2]  * **.  The 
model  neglects  elastic  and  strain  hardening  effects  and  assumes  that  the  stress 
tensor  Tij,  and  the  strain- rate  tensor,  Djj,  are  related  by  a  constitutive  relation 
^  of  the  form 

Tij  =  -p6ij  +  fi  (Dij-|divu  6ij),  (1) 

I  "4 

I  where  p  is  the  thermo^amic  pressure,  n  the  flow  resistance  coefficient,  and  u 

!  is  the  velocity  vector.  Given  the  required  relations  for  p  and  ,  (1)  may  be  in- 

coiporated  into  the  partial  differential  equations  expressing  conservation  of  mass, 
momentum,  and  energy  to  obtain  the  system  of  equations  governing  the  model. 

(See  Appendix). 

I  A  computer  code  (PICWICK)  for  solving  the  complex  system  of  equations  has 

I  been  developed  and  has  been  described  in  detail  elsewhere  [3,  4] ,  The  program 

I  is  capable  of  using  various  choices  for  the  equation  of  state  p  -  g  (p ,  I),  and  the 

'  flow-resistance  coefficients,  ju . 

;  The  principle  uncertainty  in  the  visco-plastic  model  is  the  lack  of  data  on 

(  strain-rate  effects  under  the  extreme  conditions  operative  during  hypervelocity 

impact.  It  is  necessary  to  extrapolate  data  obtained  under  much  less  severe  con- 
r  ditions  to  estimate  ji .  Part  I  of  this  report  describes  the  effect  that  these  extra- 

I  pointed  estimates  for  ju  has  on  the  structure  of  one-dimensional  stress  pulses. 

Estimates  for  fx  for  iron  and  aluminum  are  presented  which  are  compatible  both 
'  with  observed  strain-rate  effects  in  plastic  waves  and  the  very  abrupt  shock  pro¬ 

files  known  to  exist  at  pressures  greater  than  about  0.2  mb*.  These  choices  of  jx 
have  been  incorporated  into  PICWICK  and  some  calculations  made  for  axisymmetric 
impact  of  a  right  circular  cylinder  into  a  thick  target  of  like  material.  It  is  found 
that  flow  field  is  essentially  undisturbed  by  the  inclusion  of  )i>  0  in  the  early  stages 
of  the  cratering  process  for  which  calculations  with  the  full  visco-plastic  model  have 
been  completed. 

Simultaneous  with  the  investigations  of  the  strain-rate  effects,  the  PICWICK 
code  was  also  used  for  the  special  case  =  0  in  order  to  resolve  the  outstanding 
j  controversy  of  the  hypervelocity  impact  field  today.  Both  experimental  and  theo- 

I  retical  investigators  are  divided  as  to  the  manner  in  which  the  depth  of  penetration 

*  Numbers  in  square  brackets  refer  to  the  references  listed  in  the  Bibliography. 

**  The  gram-centimeter-microsecond  system  of  units  is  used  throughout  this  report, 
j  The  unit  of  pressure  becomes  the  megabar  (mb). 
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of  a  projectile  varies  with  impact  velocity  [  5] ,  In  Part  II  the  results  of  PICWICK 
calculations  for  the  Impact  of  a  cylindrical  projectile  into  a  thick  target  of  like  ma¬ 
terial  are  presented  for  lead  and  aluminum.  The  calculations  are  for  impact  veloci¬ 
ties  of  Vq  =  .76  and  2.0  cm/ /i  sec  and  are  carried  to  the  point  that  strain-rate  and 
strength  effects  become  important  (tenths  of  megabars).  It  is  found  that  momentum 
scaling  cannot  be  justified.  On  the  contrary,  it  is  shown  that  energy  scaling  is  ac¬ 
curately  predicted  if  account  is  taken  of  the  different  time  origins  when  comparing 
geometrically  similar  impact  situations  of  equal  kinetic  energy. 

In  Part  III  the  cratering  and  penetration  process  for  thin  plate  targets  is  studied. 
PICWICK  calculations  are  presented  for  iron  and  aluminum  projectiles  striking  tar¬ 
gets  of  like  material  Vq  =  .76  and  2, 0  cm/^  sec.  The  resulting  flow  field  character¬ 
istics  are  discussed  in  terms  of  the  interaction  of  the  pressure  pulse  with  the  free 
surfaces  of  the  projectile-target  system  and  with  the  rarefaction  waves  originating 
at  these  free  surfaces.  The  total  forward  momentum  of  the  projectile-target  debris 
after  penetration  is  greater  than  the  original  projectile  momentum. 


PART  I;  STRAIN-RATE  EFFECTS 


In  order  that  (1)  provides  a  constitutive  relation  for  the  entire  penetration 
process  both  p  and  ^  must  be  specified.  Analytical  expressions  for  the  thermo¬ 
dynamic  pressure  as  a  function  of  the  density,  p,  and  specific  internal  energy,  I, 

P  =  g(p.I).  (2) 

have  been  constructed.  These  functions  are  available  for  most  metals  and  cover 
the  range  of  Interest  in  blast-loading  and  hypervelocity  impact  applications  (ifp>o). 

The  bulk  of  the  studies  with  attention  to  strain-rate  effects,  however,  have 
been  made  in  which  longitudinal  plastic  stress  waves  of  relatively  low  amplitude 
are  propagated  along  bars.  The  maximum  strain-rates  obtained  are  of  order 
10"3  to  10"2  i/psec  and  the  resulting  data  have  been  fitted  into  a  polytropic  frame¬ 
work.  Several  empirical  relations  between  stress  and  strain-rate  have  been  pro¬ 
posed  which  have  proved  adequate  under  conditions  similar  to  those  in  which  the 
data  were  generated.  These  rate-dependent,  rigid-plastic  relations  may,  in  the 
absence  of  more  realistic  information,  reasonably  be  applied  in  the  later  stages 
of  flow  in  the  applications  of  interest  here. 

This  part  of  the  report  is  concerned  with  the  effect  that  the  Inclusion  of  these 
estimates  for  p  may  have  at  various  stages  of  the  flow  process.  The  important 
phases  of  the  cratering  and  penetration  process  are  illustrated  in  Figure  1, 


ONE  DIMENSIONAL  FLOW 

In  a  recent  paper  [2]  some  of  the  available  e3q)erimental  data  were  brought 
together  within  a  framework  general  enough  to  permit  the  extrapolation  of  these 
one-dimensional  empirical  formulas  in  a  natural  manner  to  the  tri-axial  case. 
Three  of  the  expressions  considered  are  as  follows: 


(Linear  Law) 


(3a) 


(Power  Law) 


(3b) 
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(3c) 


(Logarithmic  Law)  fi 


S 
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1  +  in 


provided  1t|  ^  Sq  ,  ®  in  each  case  when  |  t|  <  Sq.  Here,  and  are  the 

second  invariants  of,  respectively,  the  strain-rate  deviator  tensor  and  the  stress 
deviator  tensor.  In  (A-5)  of  the  Appendix  the  expressions  for  and  in  terms 
of  the  velocity  components  are  listed  for  the  axially  symmetric  case.  rj^,  and 
6  are  material  constants. 


In  the  case  of  a  plane  shock  propagating  in  the  z  direction,  the  relations  (A-5) 
are  simplified  by  equating  to  zero  u  and  all  derivatives  with  respect  to  u.  In  this 
part  of  the  report  in  which  only  plane  waves  are  considered  the  z  component  of 
velocity  is  denoted  by  w  rather  than  v.  Then,  =  4/3  (9w/5z)^  and  the  viscosity 
coefficients  in  (3)  reduce  to 


Approximate  values  of  the  material  parameters,  Sq,  T]t)»  6  for  the  case  of 

mild  steel  have  been  estimated  by  Perzyna  [6]  and  are  listed  in  Table  I.  Symonds 
and  Ting  [?]  have  suggested  ranges  of  values  of  Sq,  rfoi  and  6  for  the  power  law 
for  both  mild  steel  and  aluminum  alloys.  Their  estimates  are  listed  in  Table  H. 

Consider  the  plane  shock  wave  propagating  in  the  positive  direction  through 
material  of  infinite  extent,  Figure  2a,  Let  the  velocity  of  the  steady  shock  profile 
in  laboratory  coordinates  be  denoted  by  ,  and  let  w  be  the  velocity  of  the  com¬ 
pressed  material  at  any  point  in  the  profile.  Then  the  longitudinal  normal  stress 
at  a  point  on  the  profile  is  denoted  by  (‘^  the  transverse  normal  stresses 
are  Tj.j;  =  Tqq  : 


4 


(5a) 


RANKINE-HUGONIOT  DATA 


Li  a  coordinate  system  fixed  In  the  moving  profile  the  point  A,  in  the  undis¬ 
turbed  material  ahead  of  the  shock  is  fixed  relative  to  any  given  point,  B,  in  the 
shock  transition  region  (see  Figure  2a).  The  states  of  the  material  at  points  A 
and  B  may,  therefore,  be  related  by  applying  the  three  conservation  conditions: 


(mass)  (Uj- w)p=  Uj  (6) 

(momentum)  p^  -  p^  (U^-w)^  =  (7) 
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(8) 


(energy)  7  =  J  P  (Uj- w)®  +  pi  (U^- w)  (U^-  w). 

Here  the  specific  internal  energy,  I,  is  taken  to  be  zero  ahead  of  the  shock. 


From  (6)  and  (7)  one  can  deduce  that 

=  P«  w  U.  (9) 

ZZ  *^0  I 

1-X  =  w/Uj  (Xspjj/p)  (10) 

and  then  combine  with  (8)  to  write 

= -I  <“> 

I  =  I  (l-X)l  (13) 


Equations  (9) ,  (10) ,  and  (12)  are  the  familiar  Rankine-Hugonoit  relations  when  ap¬ 
plied  to  the  particular  case  where  the  point  B  is  taken  at  the  crest  of  the  shock,  so 
that  w,  I,  Til  p  then  refer  to  peak  values. 

The  Rankine-Hugonoit  relations  are  basic  to  the  experimental  studies  of  one 
dimensional  stress  pulses  induced  by  metal-metal  impact  of  two  plates,  the  one 
propelled  by  a  high  e;q}losive  [s]  .  In  these  experiments  the  petdc  values  of 
particle  velocity  w,  and  the  shock  propagation  speed  are  measured.  The  cor¬ 
responding  peak  values  of  longitudinal  stress,  722*  density  ratio,  X,  and  specific 
internal  energy  I  are  then  calculated  from  the  ttoee  equations.  This  e}q)erimentally 
determined  Hugonolt  equation  of  state  data  is  used  to  derive  Information  concerning 
the  thermodynamic  behavior  of  the  material  under  adiabatic  conditions: 


f(p.I). 


(14) 
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At  the  peak  of  the  stress  profile  Sw/3z  =  0  and  hence,  by  (5), 


P  +  ^  Sq  =  f(p,  I).  (15) 

Usually  the  metal  is  treated  as  an  inviscid  compressible  material  and  the 
term  -^4/3  Sq  does  not  appear:  g(p,  I)  =  f(p,  I)  -■^4/3  S^. 

A  function  f  ( p,  I),  defined  by  (15),  that  has  been  fitted  to  the  experimental 
data  is  of  the  form  [  9] . 


f(p,I) 


__1 _ 

<l>+<f>o 


C(aj+a2C)+0{b^+C(bi+b2C) 

+  0(c/CiC)}  . 


(16) 


where 


The  specific  internal  energy  in  the  undisturbed  material  is  taken  to  be  zero. 
Values  for  the  material  parameters  a^,  a2,  bp,  bj,  b2,  Cq,  cj,  <Pq  are  listed 
in  Table  ni  for  aluminum,  iron,  and  lead. 

In  the  following  it  will  be  convenient  to  have  a  relationship  between  the  veloc¬ 
ity  of  propagation  of  a  wave  U  and  its  peak  density  ratio  X.  To  recover  this  from 
the  equation  of  state  substitute  (11)  and  (13)  into  (16)  to  obtain 


p„U^(l-X)  =  t(p,  I) 


When  (16)  is  substituted  and  the  resulting  e}q)ression  is  solved  for  U  one  obtains 


B  -  4AC 


2A 


(17) 
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where 


A  12.,  .^,2 

A  =  4  pQ  (1-X) 


B  =  T  p 
2  ^0 


c^X  (1-X)^  +  c^X^(l-X)  -  2X^ 


b2{l-X)^  +  b^X(l-X)^  +  b^X^(l-X)-2<^^X^ 


C  =  a  X  +  a  (1-X). 

1  M 


For  X  =  1  relation  (17)  Is  indeterminate  and  must  be  replaced  by 


(18) 


SHOCK  TRANSITION  EQUATIONS 

In  this  section  the  method  of  Band  [lO]  is  extended  to  obtain  the  shape  of  the 
profile  in  the  shock  transition  region  H.  It  is  convenient  to  consider  a  coordinate 
system  moving  with  the  velocity  of  the  permanent  profile.  The  particle  velocity 
in  this  system  is  denoted  by 


q  =  Uj  -  w 


(19) 


The  conservation  equations  for  mass  and  for  momentum  for  the  one-dimensional 
case  are  (from  the  Appendix) 


at 


q 


9z 


(20) 


dz 


T 

zz 


(21) 


Since  the  profile  is  invariant  in  this  coordinate  system  the  partial  derivatives  with 
respect  to  time  vanish.  Upon  multiplying  (20)  by  q  and  adding  to  (21)  one  obtains 


9 

T"’  T 

dz  ZZ 


(22) 
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Substitute  (19)  into  (6)  to  obtain 
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At  the  crest  of  the  wave,  located  in  Figure  2a  by  the  interface  between  regions 
n  and  m,  hw/dz  =  0  and  from  (5)  and  (24)  it  follows  that 


(26) 


Similarly,  at  the  crest  of  a  wave  in  which  the  density  ratio  from  foot  to  crest  is  X 
and  which  propagates  at  velocity  U, 


S  =  Js  S  +  p  (1  -X). 
0  V  0  '^o  '  ' 


Substitution  of  (26)  and  (27)  into  (25)  yields 


ii 

az 


Oq  (Uf  -  (1  -  X) 


M  = 

az 


Pq  (Uj  -  U^)  (1  -  X) 


is 

az 


exp 

Op  (Uj  -  U^)  (1  -  X) 

-1 

.  {1  ®o  . 

. 

(27) 


(28a) 


(28b) 


(28c) 


where  U  and  X  are  related  according  to  (17). 

A  simple  numerical  scheme  has  been  used  to  calculate  the  field  variables  in 
the  shock  transition  region,  II.  X  is  treated  as  the  independent  variable,  the  range 
being  broken  into  equal  J  equal  Intervals  ^  =  (1  -  Xf)/J,  Xj  =  l-  jpX(j  =  0, 

1 . J).  From  (17),  (11),  (13),  (27) 


Xf  =  Xj  Xj.i  ^ 


1 
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and  (28)  the  corresponding  values  of  Uj,  (rzz)j>  PJ>  (d4/dz)j  maybe  computed 
in  that  order. 

From  (23)  it  follows  that  aq/az  =  Uf  SX/az,  whence  (j  =  1,  . . . ,  J): 

2  Uj  AX 

2J  -  zj-i  -  +  (aq/9z)j 

By  choosing  J  lai^e  enough  (29)  provides  the  spatial  distribution  of  the  shock  var¬ 
iables  sufficiently  accurate.  The  calculations  are  started  at  Zq  =  0  and  (aq/dz)Q  =  0. 

The  thickness  of  the  shock  profile,  AP>  is  defined  by  the  sketch  in  Figure  2b. 

It  is  easy  to  see  that 


1  -  X 


J 


AP 


where  the  subscript  i.p.  denotes  that  derivative  is  evaluated  at  the  inflection  point 
of  the  shock  profile.  Hence 


1-Xj 

'  “j  (w/sz), 


(30) 


CALCULATED  STRAIN-RATE  EFFECTS 

In  Figure  3,  the  calculated  stress  profiles  are  depicted  for  mild  steel  in  the 
case  that  the  final  compression  ratio  is  Xf  =  0. 9.  The  results  are  shown  for  three 
choices  of  the  parameters  for  each  of  the  three  forms  of  n .  The  nine  choices  of 
parameters  are  those  listed  in  Table  I.  It  is  clear  from  Figure  3a  that  the  linear 
law  spreads  the  shock  too  far  to  be  realistic  in  the  megabar  range.  The  same  is 
true  for  the  power  law  if  5  =  3.  The  profile  determined  by  the  power  law  depends 
mainly  on  the  choice  of  6,  being  relatively  Insensitive  to  changes  in  fjQ.  While  the 
logarithmic  law  provides  a  steep  front  near  the  inflection  point  it  is  extremely  slow 
in  reaching  the  asjnnptotlc  value  of  0.15  mb. 

Calculations  were  also  made  for  some  parameter  choices  in  the  case  Xf  =  0.9. 
Then  (-Tzz)f  =  Pf  0.51  mb  and  a  steady  plane  wave  should  have  an  extremely 
sharp  front.  The  thickness  of  the  profiles,  computed  by  (30),  are  listed  in  Table  I. 
The  logarithmic  law  and  the  power  law  with  6  =  5  are  compatible  with  this  require¬ 
ment. 
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On  the  basis  of  these  calculations  it  is  concluded  that  of  the  three  choices  for 
U  the  power  law  with  5  =  5  provides  a  reasonable  choice  for  the  flow  resistance 
coefficient  in  the  absence  of  more  accurate  information.  The  choices  that  are 
being  used  are: 


Aluminum: 

So  =  0.00162 

riQ  =  0.003 

6  =  8 

(31a) 

Steel: 

So  =  0.00125 

7)o  “  2.7 

6  =  6 

(31b) 

These  values  lie  within  the  range  of  values  listed  in  Table  n. 

The  choices  (31)  have  been  incorporated  into  PICWICK  and  some  calculations 
made  to  investigate: 

a)  The  strain-rate  effect  in  the  cavitation  region  behind  the  shock  front 
where  both  the  pressure  and  the  strain-rate  are  less  than  is  the  front 
itself,  and 

b)  The  strain-rate  effect  on  the  pulse  Itself  after  it  has  decreased  to  the 
order  of  tenths  of  megabars  and  below. 

At  the  time  this  report  is  being  assembled,  the  calculations  are  not  complete  but 
some  conclusions  may  be  drawn. 

Aluminum-aluminum  PICWICK  calculations  with  the  choice  (31a)  for  u  have 
been  carried  out  for  impact  velocity  Vq  =  0.76  cm/usec.  The  calculations  are 
available  to  t  =  1.081  usee.  At  this  stage  the  peak  pressure  has  dropped  below 
0.3  mb  and  only  slight  changes  in  either  (a)  or  (b)  are  discernible  when  compared 
to  an  identical  calculation  with  u  =  0.  Plots  of  the  mass  distribution  at  t  =  0.209 
and  1.081  usec  for  the  case  u>  0,  are  shown  in  Figures  5a  and  5b,  respectively. 
(This  type  of  mass  point  plot  is  described  more  fully  in  Part  H. ) 

The  fact  that  the  inclusion  of  u  does  not  alter  the  flow  process  significantly 
until  the  pressures  are  well  below  1  mb  are  used  in  the  next  part  of  this  report  as 
justification  for  continuing  the  hydrodynamic  calculations  to  relatively  low  pres¬ 
sures  .  The  limit  of  the  validity  of  the  hydrodynamic  model  probably  lies  between 
0.04  and  0.3  mb  since  Jones,  Neilson  and  Benedick  ^llj  have  observed  stress 
profiles  in  explosively  generated  plane  waves  in  aluminum  and  iron  up  to  0.04  mb. 
They  have  not,  however,  determined  the  effect  of  varying  the  load  level. 
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PART  n:  IMPACT  INTO  A  THICK  TARGET 


Immediately  upon  Impact  of  a  metal  projectile  and  metal  target  at  a  iQrperveloc- 
Ity,  pressures  of  the  order  of  megabars  are  produced  near  the  contact  interface. 
Shock  fronts  originating  at  the  contact  surface  then  propagate  upward  into  the  pro¬ 
jectile  and  downward  into  the  target  (Figure  1).  From  the  initial  Instant  a  rare¬ 
faction  wave,  originating  at  the  periphery  of  the  contact  surface,  propagates  into 
the  high  pressure  region  at  a  velocity  greater  than  the  shock  velocity.  The  rare¬ 
faction  wave  thus  catches  up  with  and  degrades  the  shock  fronts  in  the  projectile 
and  target.  The  edge  rarefaction  wave  separates  the  unrelaxed  pressure  pulse  into 
two  regions  upon  propagating  inward  to  the  axis  of  symmetry.  Rarefaction  waves 
also  occur  when  the  shock  fronts  in  projectile  anu  target  reflect  at  free  surfaces  of 
the  projectile  target  configuration. 

The  shock  dynamics  during  the  early  stages  after  impact  of  a  right  circular 
cylinder  have  been  studied  in  detail  [  12] .  It  suffices  here  to  use  the  results  of  the 
study  in  Interpreting  the  flow  fields  computed  using  the  PICWICK  code.  Calculated 
points  in  space  and  time  at  which  the  pressure  pulse,  the  rarefaction  waves,  and 
the  geometry  interact  are  listed  in  Table  IV.  The  definitions  of  the  notations  used 
in  the  following  discussion  are  also  given  there . 

After  describing  the  salient  features  of  the  flow  field,  the  total  forward  momen¬ 
tum,  the  total  outward  momentum,  and  peak  pressure  at  various  stages  of  the  flow 
process  are  studied.  It  will  be  shown  that  upon  taking  into  account  the  different 
time  origins  of  equivalent  pressure  pulses,  energy  scaling  is  accurately  predicted. 


FLOW  FIELDS  FOR  THICK  TARGETS 

In  the  basic  numerical  scheme  used  for  the  PICWICK  calculations  the  projectile 
and  target  materials  are  represented  by  discrete  mass  points .  The  distribution  of 
the  points  at  any  stage  of  the  calculations  represents  the  projectile  target  config¬ 
uration  at  that  instant.  The  area  density  of  the  points  illustrates  the  relative  mass 
density  of  the  material  and  is  indicative  of  the  pressure  existing  at  that  point  in 
space  and  time.  In  Figure  4  the  distribution  of  the  mass  point  representation  of 
the  projectile -target  configuration  at  the  Instant  of  impact  is  presented.  The 
PICWICK  code  provides  for  automatic  plots  at  the  desired  time  cycles  during  a 
computational  run.  A  sequence  of  mass  point  distributions  at  various  stages  of  a 
given  impact  calculation  illustrates  the  flow  process  in  a  vivid  manner. 
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Li  Figures  5,  6,  and  7,  representative  sequences  of  mass  particle  plots  are 
reproduced.  An  outline  of  the  original  projectile -target  configuration  is  super¬ 
imposed  on  the  plots  for  reference  purposes .  All  these  calculations  were  carried 
out  for  a  cylindrical  projectile  of  radius  equal  to  length  and  a  thick  target 
(L  =  f  =  .26192  cm  and  h  =  “  in  Figure  1).  The  volume  of  the  projectile  is  the 
same  as  that  of  a  3/32"  diameter  sphere. 

The  characteristic  mechanisms  familiar  from  experimental  studies  are  clearly 
illustrated  by  these  plots .  These  include  the  formation  of  the  backsplash  at  the 
periphery  of  the  forming  cavity  and  the  propagation  of  a  nearly  hemispherical  shock 
wave  into  the  target.  Some  of  the  basic  but  more  subtle  effects  that  have  not  been 
observed  experimentally  may  also  be  detected. 

Figure  5a  depicts  the  situation  .  209  psec  after  aluminum-aluminum  impact  at 
.  76  cm/psec.  The  rear  surface  of  the  projectile  remains  undisturbed  since  it  has 
not  yet  been  reached  by  the  shock  wave  in  the  projectile,  >  .209  according  to 
Table  IV.  The  shock  fronts  in  projectile  and  target  remain  planar  near  their  axis 
since  there  they  have  not  yet  been  caught  by  the  edge  rarefaction  wave,  tj  =  t2  >  .209. 

Figure  5b  illustrates  the  same  calculation  only  at  1.080  psec  after  impact. 

The  mesh  used  in  the  computational  scheme  has  been  repartitioned  to  permit  a 
greater  volume  of  the  target  to  be  examined  by  the  same  number  of  computational 
cells  1 3l .  At  this  stage  the  shock  in  the  target  has  been  relaxed  at  the  axis, 

1.080  >  ti,  and  propagates  on  a  nearly  hemispherical  front.  The  amplitude  of  the 
pressure  pulse,  however,  decreases  from  a  maximum  value  at  the  axis  to  zero  at 
the  surface  of  the  target  where  it  is  reflected,  causing  material  to  be  spalled  and 
ejected  backward  from  the  target  face. 

The  sequence  of  plots  in  Figure  6  again  depicts  aluminum-aluminum  impact 
but  at  Vq  =  2.0  cm/psec.  Figure  6b  depicts  a  phenomenon  that  always  occurs  in 
these  thick  target  calculations .  A  high  density  region  is  formed  in  the  target  about 
one  projectile  radius  beneath  the  center  of  impact.  From  Table  IV  it  may  be  seen 
that  for  all  cases  considered  Rj/L  =  1  and  the  high  density  region  is  located  at  the 
center  of  convergence  of  the  edge  rarefaction  wave. 

Figures  7a  through  7e  show  five  stages  of  the  flow  process  resulting  from  a 
lead  projectile  impacting  a  thick  lead  target  at  2.0  cm/psec.  Two  repartitions  of 
the  mesh  were  necessary  in  this  run. 


ENERGY  VERSUS  MOMENTUM  SCALING 

Calculations  for  like  metal  impact  of  aluminum  and  lead  are  available  (for 
Vq  =  .76  and  Vq  =  2.0  cm/psec)  from  the  instant  of  impact  until  the  pressure  has 
decreased  to  about  .  2mb.  Strain -rate  and  strength  effects  may  begin  to  become 
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important  at  these  final  pressures,  the  hydrocfynamic  model  has  been  a  valid 
approximation  to  this  point.  The  hydrodynamic  equations  may  be  put  into  a  dimen¬ 
sionless  form,  but  not  the  equation  of  state  since  it  is  an  empirical  fit.  Except  for 
this  it  would  clearly  be  possible  to  apply  the  results,  with  the  proper  scaling  factors, 
to  any  geometrically  similar  impact  situation  at  the  same  velocity  and  the  same 
material. 

To  determine  if  the  equation  of  state  form  prevents  direct  scaling  the  lead-lead 
impact  at  2.0  cm/^isec  was  re-run  with  the  configuration  scaled  up  in  dimension  by 
a  factor  of  three.  The  results  were  practically  identical  to  the  scaled  version. 

Since  the  form  of  the  equation  of  state  is  identical  for  all  the  metals  considered, 
it  was  concluded  that  scaling  is  justified  in  the  hydrodynamic  regime  for  similar 
geometries  of  a  given  metal  and  Impact  velocity. 

The  calculations  at  hand  are  sufficient  to  resolve  a  fundamental  problem  of 
iQrpervelocity  Impact  research.  Both  experimental  and  theoretical  investigators 
are  divided  as  to  the  manner  in  which  the  depth  of  penetration  of  a  projectile  varies 
with  impact  velocity  [  5) .  Some  have  concluded  that  the  cratering  process  is  de¬ 
termined  by  the  momentum  of  the  impacting  projectile;  others  have  concluded  that 
the  projectile  energy  is  the  governing  factor. 

These  hypotheses  will  be  tested  by  comparing  the  impact  calculations  at 
Vq  =  2 . 0  cm/;isec  with  the  impact  calculations  at  Vq  =  .76  cm/psec,  where  the 
latter  are  scaled  to  give  either  equal  energy  or  equal  momentum: 

(Equal  Energy) 

L  =  i  =  .26192  cm  v^  =  2.0cm/psec 

l7l  =  (Vq/vo')^^^  =  1.906  Vo'  =  .76cm/fisec  (32a) 

Mass’ /Mass  =  (L'/L)^  =  =  6.9252 

(Equal  Momentum) 

L  =  f  =  .26192  cm  Vq  =  2.0  cm/psec 

L’/L  =  (Vq/vq')^^^  =  1.38  Vq'  =  .76cm/psec  (32b) 

Mass ’/Mass  =  (L'/L)3  =  =  2.632 
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The  corresponding  formulas  for  the  depth  of  penetration  are 
(Equal  Energy) 

=  K  V  (33a) 

L  e  o 

(Equal  Momentum) 

P  1  /  o 

"r=  K  v/'*'  (33b) 

L  mo 

Two  important  characteristics  of  the  flow  field  are  the  total  forward  momentum 
and  the  total  outward  radial  momentum  present  within  the  target  at  any  instant  of 
time: 


Z+  =  E  mJ  (v^.)J  H+  =  E  (u^){  .  (34) 

i.  j  i.  j 

The  plus  indicates  that  the  sums  are  extended  only  over  positive  values  of  u  and  v. 
The  results  at  Vq  =  .76  cm/psec  under  the  hypothesis  of  equal  energy  are  scaled 
by  multiplying  R+  and  Z+  by  6.9252  and  length  and  time  by  1.906.  Momentum 
scaling  requires  that  R+  and  Z+  be  multiplied  by  2.632  and  that  length  and  time  be 
multiplied  by  1.38. 

In  Figure  8  the  calculated  total  forward  axial  momentum  is  plotted  as  a  function 
of  time  for  aluminum  impact  at  Vq  =  2.0  cm/ psec.  Also  shown  for  comparison  are 
the  total  forward  axial  momentum  when  the  vq  =  .76  cm/jitsec  results  are  energy 
scaled  and  when  they  are  momentum  scaled.  In  Figure  9  the  calculated  total  out¬ 
ward  radial  momentum  is  shown  for  the  same  three  impact  situations . 

In  Figures  10  and  11,  respectively,  the  calculated  Z+  and  R+ values  are  shown 
as  a  function  of  time  for  lead-lead  impact  at  Vq  =  2.0  cm//isec.  The  energy  scaled 
and  momentum  scaled  versions  of  the  v^  =  .76  cm/jusec  calculations  are  again  shown 
for  comparison. 

It  is  clear  from  these  results  that  momentum  scaling  cannot  be  justified,  but 
that  energy  scaling  of  the  results  of  the  Vq'  =  .76  calculations  gives  values  of  Z+ 
and  R+  that  become  closer  to  the  Vq  =  2.0  cm/psec  values  as  time  is  increased. 

If  one  only  considers  these  total  momentum  curves,  it  appears  that  energy  scaling 
slightly  overestimates  the  velocity  dependence  of  the  hypervelocity  damage  mech¬ 
anism  and  that  the  exponent  in  the  penetration  formula,  (33),  should  be  slightly 
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smaller  than  2/3.  This  was  the  conclusion  of  Walsh  and  Tillotson  [is]  in  their 
recent  and  important  paper.  They  found  that  an  exponent  between  .  58  and  .  72 
brought  their  energy  equivalent  Z+  and  R+  curves  into  close  agreement.  Different 
values  of  the  exponent  were  required  for  different  velocity  ranges  and  for  different 
metals. 

The  results  discussed  in  the  next  section  show  that  this  apparent  variation  in 
the  exponent  for  different  materials  and  velocity  ranges  may  be  e^qplained  by  the 
failure  to  account  for  a  time  differential  when  comparing  equivalent  energy  impacts. 


THRESHOLD  OF  ENERGY  SCALING 

When  comparing  the  deformation  waves  associated  with  two  impact  velocities 
of  equivalent  energy,  v^  >  v,  the  characteristics  during  the  initial  stages  are  of 
course  quite  distinct.  The  higher  velocity  projectile  will  produce  larger  pressures 
upon  impact  and  the  pressure  pulse  will  propagate  at  a  faster  velocity  into  the  tar¬ 
get.  The  pressure  pulse  associated  with  the  lower  velocity  will  require  more  time 
to  reach  a  specified  distance  below  the  target  surface. 

To  determine  the  time  lag,  1 1,  between  pressure  pulses  of  equal  energy  the 
peak  axial  pressure  has  been  studied.  In  Figure  12  the  distance  of  the  peak  axial 
pressure  below  the  original  target  surface  is  plotted  as  a  function  of  the  time  for 
aluminum-aluminum  impact  at  Vq  =  2. 0  cm/jU  sec  and  for  the  energy  equivalent 
situation  for  impact  at  Vq  '=  .  76  cm/ju  sec.  In  Figure  13  the  corresponding  peak 
axial  pressure  versus  time  curves  are  shown  for  lead-lead  impact.  The  time  lag 
for  the  aluminum  and  the  lead  are  ^t  =  .  35  and  At  =  .  75  ^  sec,  respectively. 

Since  the  pressure  pulse  is  the  basic  cause  of  the  cratering  process  two  im¬ 
pact  situations  should  be  compared  at  times  when  the  pulses  are  equal  distance 
below  the  target  surface.  Accordingly,  the  energy  scaled  versions  of  Z+  and  R+ 
in  Figures  8  and  9  should  be  shifted  to  the  right  a  distance  corresponding  to 
At  =  ,  35  fisec,*  The  corresponding  curves  in  Figures  10  and  11  are  to  be  shifted 
At  =  .75  jusec  to  the  right.  In  each  figure  the  dashed  curve  represents  the  energy 
scaled  version  with  the  appropriate  time  correction.  For  both  aluminum  and  lead 
this  correction  leads  to  excellent  agreement  between  the  energy  equivalent  Impact 
situations. 

In  Figures  14  and  15  the  magnitude  of  the  peak  axial  pressure  pulse  as  a  func¬ 
tion  of  the  distance  below  the  target  surface  is  shown  for  Vq  =  2. 0  cm/ n  sec  and  for 
the  momentum  and  energy  equivalent  versions  of  the  Vg '  =  .76  cm/jU  sec  calculations. 
Figure  14  illustrates  the  aluminum  calculations  and  Figure  15  shows  the  lead  calcu¬ 
lations.  For  aluminum  impact  energy  scaling  results  in  overlapping  values  for  the 
pressures  when  p  <  .2  mb.  The  results  are  especially  convincing  for  the  lead  cal¬ 
culations  where  the  energy  scaling  gives  values  that  coalesce  for  p  <  1. 0  mb. 

♦The  same  shifting  is  required  for  the  momentum  scaled  version  but  would  only  in¬ 
crease  the  disparity  between  the  momentum  equivalent  situations. 


One  megabar  Is  &r  above  the  threshold  of  strength  and  strain-rate  effects  in 
lead  and  hence,  penetrations  must  certainly  follow  the  two -thirds  power  law,  (32a). 
On  the  other  hand  it  is  uncertain  that  0.2  mb  is  greater  than  the  threshold  pressure 
for  strength  and  strain-ratio  effects  for  aluminum  or  aluminum  alloys.  However, 
if  experimental  data  shows  that  the  penetration  of  geometrically  similar  projectiles 
into  their  targets  of  a  given  alloy  varies  according  to  Vq  2/3  at  the  upper  end  of  the 
experimental  range,  it  is  safe  to  extrapolate  at  least  to  Vq  =  2.0  cm/ fi  sec. 

A  corollary  to  this  is  that  if  the  2/3  law  is  attained  for  two  aluminum  alloys 
but  with  different  values  of  Ke,  this  difference  will  persist.  Strength  and  strain- 
rate  effects  will  not  disappear  with  increasing  impact  velocity.  They  will  continue 
to  become  predominant  during  the  latter  stages  of  the  cratering  process. 
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PART  ni;  IMPACT  INTO  PLATES 


Since  practical  structures  are  not  of  unlimited  thickness,  impacts  into  plate 
targets  have  also  been  Investigated.  The  presence  of  the  opposite  free  surface  of 
the  target  allows  the  relief  of  the  pressure  pulse  by  a  reflected  rarefaction  wave 
resulting  in  the  formation  of  the  familiar  spallation  bubble,  Figure  1.  Relatively 
thin  plates  (thickness  less  than  half  the  ballistic  limit)  are  completely  penetrated 
before  the  pressure  pulse  has  attenuated  sufficiently  for  strength  and  strain-rate 
effects  to  become  important.  These  will  determine  the  final  diameter  of  the  hole 
in  the  target,  but  a  description  of  the  debris  in  the  spallation  bubble  may  be  ob¬ 
tained  from  the  hydrodynamic  model. 

The  calculations  were  all  made  for  a  projectile-target  configuration  with 
h  m  t  aha  .26192  cm,  but  as  discussed  earlier  the  results  may  be  scaled  to 
any  geometrically  similar  impact  of  the  same  materials  at  the  same  velocity. 


FLOW  FIELDS  FOR  PLATE  TARGETS 

The  sequence  of  Figures  16a  through  16d  depict  the  pressure  pulse  propagation, 
the  cratering  process  and  the  creation  of  a  spallation  bubble  for  aluminum-aluminum 
Impact  at  Vq  *  .76  cm/ju  sec.  Figures  17a  and  17b  depict  aluminum-aluminum  im¬ 
pact  at  2.0  cm/iisec  and  Figures  18a  through  18d  show  the  results  of  the  calcula¬ 
tions  for  iron-iron  impact  at  .76  cm/|isec.  An  outline  of  the  projectile-target 
configuration  at  the  instant  of  impact  is  superposed  on  each  plot  to  provide  a  frame 
of  reference. 

The  pressure  and  velocity  fields  are  also  depicted  for  the  Vq  ■  .76  cm//isec 
calculations  for  both  the  aluminum-aluminum  impact  (Figures  19a  through  19e)  and 
the  iron-iron  impact  (Figures  20a  through  20e.) 

When  a  shock  wave  reaches  a  free  surface  it  must  be  reflected  in  such  a  man¬ 
ner  that  the  stress  normal  to  the  surface  is  zero.  For  normal  incidence  this  is 
accomplished  by  the  generation  of  a  reflected  wave  of  amplitude  equal  to  the  inci¬ 
dent  wave,  but  of  opposite  sign;  and  by  accelerating  the  reflected  surface  to  twice 
the  particle  velocity  of  the  incident  wave.  By  symmetry  the  shock  waves  in  the 
projectile  and  target  must  reach,  respectively,  the  rear  of  the  projectile  and  the 
opposite  surface  of  the  target  at  normal  incidence  along  the  axis.  Moreover,  from 
Table  IV,  it  may  be  seen  that  for  the  aluminum-aluminum  impacts  both  at  Vg  »  .76 
and  2.0  cm/^sec,  I  <  Iq  and  h  <  hg;  consequently  there  will  be  an  axial  region 
in  which  the  shock  fronts  have  not  been  relaxed  by  the  edge  rarefaction  wave.  For 
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the  iron-iron  impact  the  shocks  at  the  axis  will  be  caught  by  the  edge  rarefaction 
just  prior  to  the  time  the  shocks  reach  the  free  surface. 

At  the  opposite  surface  of  the  target,  Initially  at  rest,  the  particle  velocity  be¬ 
hind  the  shock  is  Vq/2  so  that  after  the  pressure  pulse  reaches  the  surface  the 
particle  velocity  becomes  Vq.  This  sudden  increase  in  particle  velocity  is  illus¬ 
trated  by  comparing  Figure  19b  with  Figure  19c  and  Figure  20b  with  Figure  20c. 
For  points  away  from  the  axial  region  the  shock  will  have  been  attenuated  by  the 
edge  rarefaction  prior  to  reaching  the  free  surface.  The  particle  velocity  will 
be  Vq/2  -v(r)  where  the  velocity  lag,  v{r)  >  0,  increases  with  radial  distance  r. 
After  the  reflection  of  the  pulse  at  r,  the  particle  velocity  is  Vo-2v(r)  so  that  the 
velocity  lag  is  doubled. 

At  the  rear  surface  of  the  projectile  the  particle  velocity  in  the  central  region 
is  -Vo/2  in  a  coordinate  system  fixed  in  the  rear  surface.  After  reflection  the 
particle  velocity  in  this  system  is  -v©.  In  the  laboratory  system,  however,  the 
new  particle  velocity  is  zero.  This  sudden  decrease  in  the  particle  velocity  is 
Illustrated  by  comparing  Figure  19b  with  Figure  19c  and  Figure  20b  with  Figure 
20c.  It  is  also  demonstrated  by  comparing  Figure  16b  with  Figure  16c  and  Figure 
18b  with  Figure  18c.  The  latter  plots  show  that  the  rear  mass  particles  in  the 
projectile  have  remained  nearly  stationary.  At  points  away  from  the  axial  region 
the  particle  velocity  after  the  reflection  becomes  2v(r),  l.e.,  increases  with  r. 

The  axial  portion  of  the  spallation  bubble  is  increasing  in  length  at  the  rate  Vq 
while  the  outer  portion  is  increasing  in  length  at  the  rate  Vq  -  4v(r) .  This  produces 
the  elongated  geometry  of  the  bubble  that  is  observed  experimentally. 


RELEASE  OF  PRESSURE  PULSE 

One  inference  of  the  calculations  is  that  for  the  plate  thickness  and  projectile 
velocities  studied,  the  penetration  is  completed  before  much  of  the  energy  carried 
by  the  pressure  pulse  can  be  transmitted  radially  into  the  target.  This  is  demon¬ 
strated  by  the  pressure  profiles  shown  in  Figures  19a  through  20e.  The  pressure 
in  the  central  region  of  the  projectile-target  configuration  decreases  to  zero  soon 
after  the  pressure  pulse  is  relieved  by  the  presence  of  the  opposite  free  surface  of 
the  target.  Once  the  pressure  in  the  central  region  decreases  to  zero  the  pulse  that 
has  propagated  radially  into  the  target  is  no  longer  supported  and  attenuates  rapidly. 
For  both  aluminum  and  iron  impact  at  Vq  =  .76  cm/ju  sec  the  maxinum 
pressure  has  decreased  to  approximately  .3mbatr  =  2L.  Strength  and  strain- 
rate  effects  will  play  an  important  role  in  the  subsequent  radial  deformation  of  the 
target. 

The  sudden  release  of  the  pressure  pulse  results  in  a  rapid  expansion  of  the 
material  accompanied  by  a  decrease  in  the  internal  energy  and  a  corresponding 
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increase  in  the  kinetic  energy.  This  dramatic  change  in  the  partition  of  the  energy 
is  depicted  in  Figures  21a,  22a,  and  23a  for,  respectively,  aluminum  Impact  at 
.76  and  2.0  cm/^tsec  and  iron  Impact  at  .76  cm//isec.  Curves  are  also  shown  de¬ 
picting  the  partition  of  energy  for  the  corresponding  thick  target  impact  situations. 
Even  after  the  sudden  decrease,  however,  the  total  internal  energy  reaches  a 
plateau  representing  the  release  temperature  resulting  from  the  increase  in  en¬ 
tropy  produced  by  the  shock. 

In  the  special  case  of  two  semi-infinite  bodies  impacting  at  velocity  Vq,  the 
kinetic  energy  and  internal  energy  behind  the  shock  front  are  equal,  I  =  k.e.  * 
Vq^/S.  This  equlpartitlon  represents  an  upper  limit  on  the  value  for  internal  energy 
since  there  are  no  rarefaction  waves  present  in  such  an  idealized  situation. 

Figures  21,  22,  and  23  also  show  the  total  forward  momentum  and  the  total 
outward  radial  momentum  for  all  three  impact  situations.  Both  are  seem  to  increase 
with  time  but  at  a  much  slower  rate  than  the  corresponding  thick  target  impact 
once  the  pressure  pulse  has  been  relieved. 

Since  and  R).  are  important  parameters  when  considering  the  effectiveness 
of  plate  targets  as  shields  for  a  protected  structure,  it  may  be  concluded  that  the 
plate  should  only  be  as  thick  as  required  to  break  up  the  projectile.  If  the  spacing 
is  large  enough  the  expanding  bubble  will  then  be  dispersed  over  a  wide  area  of  the 
structure. 
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CONCLUSIONS 


The  most  important  conclusion  of  this  report  is  that  the  craters  produced  by  the 
impact  of  geometrically  similar  projectiles  of  the  same  metal  into  a  thick  target  of 
like  metal  will  be  the  same  if: 

(1)  The  projectiles  are  of  equal  kinetic  energy 

(ii)  The  impact  velocity  exceeds  a  threshold  value  for  energy  scaling 

This  has  been  concluded  from  a  detailed  study  of  aluminum-aluminum  and  lead- 
lead  Impact  at  Vq  »  0.76  and  2.0  cm/jisec.  The  calculations  were  carried  out 
for  cylindrical  projectiles  of  radius  equal  to  length,  L  *  -t . 

The  formula  for  the  depth  of  penetration  corresponding  to  (i)  and  (li)  is 


(Vq  ^  V  *), 


where  Vq*  is  the  threshold  velocity  for  the  2/3  law  for  a  particular  metal  and  Kg  the 
corresponding  experimental  constant.  These  parameters  are  known  from  experi¬ 
ment  to  vary  among  alloys  of  nominally  the  same  metal,  reflecting  strength  and 
strain-rate  effects.  The  calculations  indicate  Vq  <  .76  for  lead  and  for  aluminum 
except  possible  aluminum  alloys  in  which  strain-rate  and  strength  effects  persist 
for  pressures  up  to  .2  mb. 

These  conclusions  are  based  on  a  very  strict  criterion  for  equivalence  of  two 
impact  cratering  processes.  Since  the  pressure  pulse  in  the  target  is  the  basic 
cause  of  the  cratering  mechanism,  the  two  processes  are  compared  at  times  when 
the  locations  of  their  peak  pressures  are  an  equal  distance  below  the  target  sur¬ 
face.  Based  on  the  time  since  impact,  the  process  associated  with  the  slower  im- 
,pact  velocity  Vg '  lags  the  other,  Vq  ,  by  At.  For  equivalence  it  is  required  that  at 
these  corresponding  times  the  total  forward  axial  momentum  Z+,  the  total  outward 
radial  momentum,  R^,  and  the  amplitude  of  the  pressure  pulse  become  equal  prior 
to  the  onset  of  strength  and  strain-rate  effects. 

If  one  simply  compares  Z+  and  R.).  on  the  basis  of  time  since  impact,  ignoring 
the  time  lag  At,  an  exponent  somewhat  less  than  n='2/3  would  be  estimated  for  the 
velocity  dependence  of  the  penetration  law.  Since  At  increases  with  increasing 
Vo -Vo' ,  the  apparent  value  for  the  exponent  would  decrease  with  increasing  impact 


velocity.  This  is  believed  to  explain  the  discrepancy  between  the  present  conclu¬ 
sions  and  those  of  Walsh  and  Tillotson  [13]  who  estimated  n  ■  .  65  (.  5  <  Vq  <  1), 
n  ■  .62  (1  <  Vq  <  4),  andn  »  .61  (10  <  Vq  <  25). 

Calculations  were  also  performed  for  aluminum  and  iron  plate  targets  im¬ 
pacted  by  like  metal  projectiles  at  Vq  ■  .76  and  2.0  cm/jusec  (plate  thickness  h, 
h^L^-t).  The  pressure  pulse  is  relieved  by  the  spallation  at  opposite  free  surface 
before  much  of  the  energy  can  be  transmitted  radially  into  the  target.  The  pres¬ 
sure  in  the  central  part  of  the  configuration  rapidly  falls  to  zero  and  an  expanding 
spallation  bubble  is  formed.  The  material  in  the  bubble  is  completely  shattered 
and  will  remain  at  a  high  temperature.  The  total  forward  momentum  of  the  con¬ 
figuration  after  impact  exceeds  the  original  axial  momentum.  If  the  thickness  of 
the  plate  is  less  than  the  projectile  radius  h<  L,  the  leading  central  portion  of  the 
spallation  bubble  travels  at  the  velocity  of  the  projectile  at  impact. 

These  results  indicate  that  the  use  of  thin  plates  as  shields  for  structures  could 
result  in  greater  damage  to  the  structure  if  the  spacing  between  shield  and  struc¬ 
ture  is  not  great  enough.  The  meteor  bumper  problem  will  be  studied  during  the 
next  six  months  of  the  contract. 
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APPENDIX 


Let  p,  p,  u  =  (u,  0,  V),  and  I  denote  the  density,  pressure,  velocity,  and 
specific  internal  energy,  respectively.  Then  the  axisymmetric  formulation  of 
the  visco-plastic  equations  may  be  written  in  the  form: 

(Mass)  |^+  v-^+  pdivu=0  (A-1) 


.T,  j,  1 ,,  .  .  ^  ^  L  r,  \  ^rr  ^66  ^rz  ,, 

(Radial  Momentum)  P^+u-^+v-^J  =  -^  iP+sJ  +  — - —  +  —  (A-2) 


(Axial  Momentum)  ±  j  (A-3) 


'I 


/dl  bl  dl\  -  9 

(Energy)  ^  “ir  ''W  P  • 


(A-4) 


Here  the  flow  is  assumed  to  be  strictly  adiabatic  and  the  following  notations  have 
been  introduced: 


P  =  -p  p  div  u 
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TABLE  I 


Values  for  mild  steel  of  the  material  parameters  occurring  in  relations  (3).  These 
values  were  estimated  by  Perzyna  [6] .  Also  listed  are  the  corresponding  shock 
thicknesses  calculated  for  the  indicated  final  compression  ratios.  The  units  are  in 
the  gm-cm-^sec  system. 


Form 

No. 

S  X  10^ 

0 

V 

0 

6 

<1 

of 

M 

(Xj=.9) 

n 

o 

00 

Linear 

B 

1.62 

5.2 

B 

10.3 

Law 

2 

1.62 

3.9 

1 

13.0 

3 

1.62 

3.1 

B 

17.4 

7.32 

Power 

4 

1.62 

23.1 

5 

.01 

3.0 X  lO"® 

Law 

.0 

5 

1.62 

5.2 

3 

.20 

2.2x  10 

6 

1.62 

3.9 

3 

.15 

— 

Logarithmic 

B 

1.62 

11.7 

B 

.034 

__ 

Law 

8 

1.62 

7.8 

■ 

.022 

~ 

1.62 

1 

.011 

3.6  X  lo"^^ 
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TABLE  n 


Ranges  of  values  estimated  Iqr  Symonds  and  Ting  [7]  for  the  parameters  in  the 
power  law  for  mild  steel  and  structural  aluminum.  The  estimates  are  based  on 
the  results  of  a  number  of  experimenters. 


Metal 

S  X  10^ 

0 

^0 

6 

Mild  Steel 

1.25 

2.7  -  18 

CO 

1 

A1  6061- T6 

1.62 

. 003  -  .  154 

4-8 

TABLE  m 


Equation  of  state  parameters,  occurring  in  equation  (16),  expressed  in  the  gm-cm-jusec 
system  of  units. 


Metal 

"i 

"2 

D 

*^2 

<t> 

^0 

B 

A1 

1.1867 

.763 

- 1 

3.4448 

1.5451 

.9643 

.43382 

.  54873 

m 

2.702 

Fe 

7.78 

31.18 

9.591 

15.676 

4.634 

.3984 

.5306 

9.0 

7.86 

g 

1.4844 

1. 6765 

8.7317 

.96473 

2.6695 

.27732 

.43709 

3.3 

11.34 
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TABLE  IV 


Points  in  time  and  space  at  which  the  shock  fronts,  rarefaction  waves,  and  the  con¬ 
figuration  surfaces  Interact.  ThevalueslnthistableweresupplledbyDr.  J.  F.  Heyda 
and  are  valid  for  like  metal  impact. 


Point  of 

V  = 
0 

.76  cm/fi 

sec 

V  = 
0 

2.0  cm/^sec 

hiteraction 

A1 

Fe 

Pb 

A1 

Fe 

Pb 

tj/L  =  t^/L 

1.116 

.963 

1.208 

.658 

.499 

.737 

Rj/L 

1.178 

.963 

.938 

1.188 

.945 

1.142 

R/L 

.330 

.231 

.019 

-.127 

-.053 

-.332 

tj/L 

.891 

.827 

1.090 

.581 

.456 

.683 

R/L 

.339 

.314 

.414 

.581 

.456 

.683 

TjA 

.947 

1,000 

1.289 

.553 

.528 

.645 

^  /L  =  h  /L 

0  0 

1.178 

.963 

.938 

1.188 

.945 

1.142 

t  j  =  time  for  edge  rarefaction  to  catch  shock  front  in  target  at  the  axis 
t„  =  time  for  edge  rarefaction  to  catch  shock  front  in  projectile  at  the  axis 
Rj  =  position  of  shock  front  in  target  at  t^^ 

R  =  position  of  shock  front  in  projectile  at  t„ 

£t  u 

t^  =  time  for  edge  rarefaction  to  first  reach  axis  at  moving  interface 
R^  =  position  of  interface  at  t^ 

T  j  =  time  required  for  shock  in  projectile  to  reach  its  rear  surface 
I  =  projectile  length  for  which  t  =  t  . 

0  ^  X 

h  =  plate  thickness  for  which  edge  rarefaction  reaches  axis  in  target 
simultaneously  with  the  shock 
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Figure  3.  Calculated  shock  profiles  in  mild  steel  for  a  final  compression 
ratio  Xf  ■  0.9  when  ju  is  assumed  to  follow  (a)  linear  law, 

(b)  power  law,  (c)  logarithmic  law. 


Figure  5a.  The  mass  point  plot  .209  psec  after  impact  of  an  aluminum  projectile  into  a  thick 
aluminum  target  at  .76  cm/usec. 


aluminum  target  at  .  76  cm/ nsec . 


Figure  6a.  The  mass  point  plot  .249  ju  sec  after  impact  of  an  aluminum  projectile 
into  a  thick  aluminum  target  at  2.0  cm/^sec. 


Figure  6b.  The  mass  point  plot  .  550  sec  after  impact  of  an  aluminum  projectile 
into  a  thick  aluminum  target  at  2. 0  cm//isec. 


Figure  7a.  The  mass  point  plot  at  t  =  .  249  fiaec  after  impact  of  a  cylindrical 
lead  projectile  into  a  lead  target  at  Vq  =  2. 0  cm//xsec. 


Figure  7b.  The  mass  point  plat  at  t  =  .  525  fisec  after  impact  of  a  cylindrical 
lead  projectile  into  a  thick  lead  target  at  Vq  =  2. 0  cm/jisec. 


Figure  7c.  The  mass  point  plot  at  t  =  1.297  /xsec  after  impact  of  a  cylindrical 
lead  projectile  into  a  thick  lead  target  at  Vo  =  2. 0  cm/u  sec. 


Figure  7d.  The  mass  point  plot  at  t  =  3.086  sec  after  impact  of  a  cylindrical 
lead  projectile  into  a  thick  lead  target  at  =  2. 0  cm/jxsec. 


TOTAL  FORWARD  AXIAL  MOMENTUM  (GM~CM//aSEC) 


Figure  8.  Total  forward  momentum  in  the  target  as  a  function  of  the  time 
since  impact  of  an  aluminum  projectile  and  a  thick  aluminum 
target. 
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Figure  9.  Total  positive  radial  momentum  in  a  thick  aluminum  target  at  indi' 
cated  time  since  impact  by  an  aluminum  projectile. 
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Figure  10.  Tot 
sin 


TIME  AFTER  IMPACT  (/xSEC) 


forward  momentum  in  a  thick  target  as  a  function  of  the  time 
impact  of  a  lead  projectile  and  target. 


0  12  3  4  5  6 

TIME  AFTER  IMPACT  {ft  SEC) 


Figure  11.  Total  positive  radial  momentum  in  a  thick  target  at  indicated 
time  since  lead-lead  impact. 
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DISTANCE  FROM  TARGET  SURFACE  OF  PEAK  AXIAL  PRESSURE  (CM) 


Figure  12.  The  position  of  the  peak  axial  pressure  as  a  function  of  time  for 
.aluminum-aluminum  impact. 
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Figure  13 .  The  position  of  the  peak  axial  pressure  as  a  function  of  time  for 
lead-lead  impact. 


48 


PEAK  AXIAL  PRESSURE  (Mb) 


Figure  14.  Peak  axial  pressure  in  aluminum  shown  as  a  function  of  distance  below 
original  surface  of  a  thick  target. 
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DISTANCE  FROM  TARGET  SURFACE  (CM) 


Figure  15.  Peak  axial  pressure  in  lead  shown  or  function  of  distance  below  the 
original  surface  of  a  thick  target. 
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Projectile-target  material  distribution  ,4 
an  aluminum  projectile  into  an  aluminum 


Figure  16c.  Projectile-taxget  material  distribution  .550  ^sec  after  impact  of 
an  aluminum  projectile  into  an  aluminum  plate  at  .76  cra/fisec. 
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Figure  17b.  Projectile-target  material  distribution  at  ,410  ^sec  after  impact  of 
an  aluminum  projectile  into  an  aluminum  plate  at  2.0  cm/jxsec. 


Figure  18a.  Projectile-target  material  distribution  .255  fxsec  after  impact  of 
an  iron  projectile  into  an  iron  plate  at  .  76  cm/jisec. 


Projectile-target  material  distribution 
an  iron  projectile  into  an  iron  plate  at 


Figure  18d.  Projectile-target  material  distribution  .676  jusec  after  impact  of 
an  iron  projectile  into  an  iron  plate  at  .76  cm/fisec. 


Figure  19a.  Pressure  and  velocity  fields  at  .136  fisec  after  impact  of  an  alumi¬ 
num  projectile  into  an  aluminum  plate  at  .76  cm//isec. 


Figure  19e.  Pressure  and  velocity  fields  at  .674  fisec  afti 
num  projectile  into  an  aluminum  plate  at  .76 
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Figure  20a.  Pressure  and  velocity  fields  at  .  136  fjisec  after  impact  of  an 
iron  projectile  into  an  iron  plate  at  .76  cm/jzsec. 


Figure  20b.  Pressure  and  velocity  fields  at  .203  fisec  after  impact  of  an  iron 
projectile  into  an  iron  plate  at  .76  cm/ fisec. 


Figure  20c.  Pressure  and  velocity  fields  at  .421  ^sec  after  impact  of  an  i 
projectile  into  an  iron  plate  at  .76  cm/usec. 


igure  20(1.  Pressure  and  velocity  fields  at  .539  fisec  after  impact  of 
projectile  into  an  iron  plate  at  .76  cm/usec. 


Figure  20e.  Pressure  and  velocity  fields  at  .674  fxsec  after  impact  of  an  iron 
projectile  into  an  iron  plate  at  .76  cm/^sec. 
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Comparison  of  flow  field  characteristics  for  aliuninum-aluminum  impact  in  thick  targets 
and  in  plate  targets  at  .  76  cm//i  sec:  (a)  Energy  partition,  (b)  Total  forward  axial 
momentum,  and  (c)  Total  outward  radial  momentum. 
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SUMMARY 

The  principle  uncertainty  in  the  visco-plastic  model  of  hyper¬ 
velocity  impact  is  the  lack  of  data  on  strain-rate  effects  imder  thej 
severe  conditions  of  interest.  Estimates  of  n  for  iron  and  alumi¬ 
num  are  first  determined  that  are  compatible  both  with  observed 
strain-rate  effects  in  plastic  waves  and  the  very  abrupt  shock 
profiles  known  to  exist  at  pressures  greater  than  about .  2  megabafr. 

In  the  second  part  of  the  report  detailed  calculations  are  pre¬ 
sented  describing  the  cratering  process  for  a  cylindrical  projectil|e 
impacting  a  thick  target  of  like  metal.  It  is  shown  that  momentum 
scaling  cannot  be  justified;  energy  scaling  for  geometrically  simij 
lar  impact  situations  is  accurately  predicted  upon  taking  into 
account  the  different  rates  of  shock  propagation  during  the  early 
stages  of  the  cratering  process.  The  depth  of  penetration  is  there¬ 
fore  related  to  impact  velocity  v^  according  to  I 

Vq  =  Vq*,  Strength  effects  of  the  target  do  not  vanish  with  increasing 
impact  velocity. 

In  the  third  part  of  the  report  the  cratering  and  the  penetration 
process  for  thin  plates  is  studied.  PICWICK  calculations  are  pre-j 
seated  and  the  results  discussed  in  terms  of  the  interactions  of 
shock  waves,  rarefaction  waves,  and  configuration  geometry. 
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